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We demonstrate that the different magnetic states of a Mni2 molecule can be distinguished in a 
two-probe transport experiment from a complete knowledge of the current-voltage curve. Our re- 
sults, obtained with state-of-the-art non-equilibrium transport methods combined with density func- 
tional theory, demonstrate that spin configuration-specific negative differential resistances (NDRs) 
appear in the I-V curves. These originate from the interplay between electron localization and the 
re-hybridization of the molecular levels in an external electric field and allow the detection of the 
molecule's spin-state. 

PACS numbers: 



The possibility of addressing and manipulating single 
spins represents the ultimate limit for the magnetic data 
storage industry. It may also lead to a technological plat- 
form for the integration of logic and data storage func- 
tionalities on the same device. A crucial aspect of any 
single spin-device is the development of a protocol for 
the electrical readout of the magnetic state of the device. 
Here we explore the possibility of inferring the spin con- 
figuration of a magnetic molecule by a detailed knowledge 
of the I-V characteristic of a two-probe device having the 
molecule as resistive element. 

In particular we consider the case of Mni 2 sandwiched 
between non-magnetic electrodes and we compare the 
electrical response of two possible spin configurations of 
the molecule. First we look at the ground state (GS), 
where the spins of the eight peripheral Mn 3+ ions align 
antiparallel to those of the inner Mn 4+ to give an over- 
all S = 10 spin state [1 j. The second is a spin- flip (SF) 
configuration, in which the spin directions of one Mn 3+ 
and one Mn 4+ are reversed with respect to those in the 
GS. The total spin-projection for such a state is S = 9. 
The most notable difference between the I-V curves ob- 
tained for the different spin states is the much higher low- 
bias current of the GS configuration and the presence of 
NDRs, which are specific to the spin state. These are a 
consequence of orbital re- hybridization under bias, which 
causes a highly non-linear bias-dependent coupling of the 
molecular levels to the electrodes. We predict that this is 
a general feature of molecular junctions characterized by 
closely spaced orbital multiplets, such as those appear- 
ing in magnetic molecules. Importantly, since both the 
orbital symmetry and their localization depend on the 
molecule's spin-state, we expect different I-V fingerprints 
for different magnetic configurations. This means that 
the overall molecular magnetic configuration is readable 
entirely from a single non-spin-polarized current readout. 

Two experiments have recently measured the I-V 
curve of three-terminal devices using Mni2 as the resis- 
tive element [2 , 3 , in both cases revealing features asso- 
ciated to the magnetic state of the molecule. NDRs have 
been ascribed to selection-rule forbidden transitions be- 
tween different charging states [4j [5], with the theory 
being based on model Hamiltonians and parameters ex- 



tracted from Mni2 single crystal properties. A direct 
quantitative analysis of the results is complicated by the 
large variation in I-Vs from device to device and by the 
low device production yield. This stems from the gen- 
eral fragility of Mni2 in single molecule form on surfaces 
against fragmentation [6 and redox [7 . The spectro- 
scopical detection of Mn 2+ in a Mni 2 monolayer [7] is 
proof of such fragility. Furthermore, when the molecule 
remains intact without changing its oxidation state, still 
the magnetic response is rather sensitive to the surface 
environment and different from that of a single crystal. 
Most notably the anisotropy disappears [8], making the 
interpretation of the transport features in terms of selec- 
tion rule-forbidden transitions more complicated. Never- 
theless, these pioneering transport measurements provide 
an indirect evidence for the different spin states of the 
molecule. Interestingly such evidence is obtained without 
using spin-polarized currents as in spin-polarized scan- 
ning tunnel microscopy (STM) [9j [10] and so may repre- 
sent an intriguing alternative for single-spin detection. 

We begin our analysis by presenting the elec- 
tronic structure of the archetypal Mni2 acetate cluster 
[Mni 2 Oi2(CH 3 COO)i6(H 2 0)4] both in the GS and in 
the SF configuration. The calculations are performed 
with density functional theory (DFT) using the general- 
ized gradient approximation (GGA) [11] as numerically 
implemented in the pseudo-potential, local orbital basis 
set, siesta code [3]. Details of the numerical method 
are presented as supplementary material. The density 
of states (DOS) and the associated local DOS for a few 
of the molecular orbitals around the Fermi energy (E-p) 
are presented in figure [6j Clearly the DOS is entirely 
dominated by the Mn 3+ d shells (see also supplementary 
materials), with an energy gap between the highest oc- 
cupied molecular orbital (HOMO) and the lowest unoc- 
cupied molecular orbital (LUMO) of about 0.3 eV. This 
is considerably smaller than the estimated experimental 
gap (1.4 eV [13 ) as expected from the self-interaction 
contributions [14] to the GGA functional. However such 
a reduced gap has little qualitative relevance for the low- 
bias I-V, since the nature of the HOMO manifold is in- 
sensitive to the choice of the functional [T3] . 

In the GS the DOS is entirely spin-polarized around 
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E F with a triply degenerate HOMO followed by a singly 
degenerate HOMO-1. Such a degeneracy is lifted in 
the SF configuration where the HOMO is formed by 
four closely spaced Mn 3+ levels, one of which has the 
opposite spin. Overall the total DOS (not considering 
the spin polarization) is rather similar in the two cases, 
suggesting that an electrical measurement, which is not 
sensitive to the spin direction, will hardly be able to 
distinguish between them. However there is an impor- 
tant difference between the two configurations. This is 
the degree of localization of the various HOMO levels. 
In the GS the local DOS shows amplitude uniformly 
distributed around all the eight Mn 3+ ions, while in 
the SF state the minority HOMO is extremely localized 
around the flipped Mn 3+ ion. This has a profound 
influence over the electrical response of the molecule. 



150 

100 

■3 50 - 
a 

U 



g-50 



-100 - 
-150 



(a) 



GS1 



■ ' 1 1 1 ' ' 1 1 1 1 



Total 
Ma 



Hip 



■(b) 



SF1 



-- / 



SF2 



1 "1 1 1 "1 1 
Majority- 



Minority_ 
i i i ' I i i 1 1 1 i 



-0.4 -0,2 0.2 0.4 



-0.4 -0.2 0.2 0.4 



< < c -?C J. 

GS1 V* n SF1 




SF2 




FIG. 1: Density of states around the Fermi level (purple ver- 
tical line) for [Mni 2 Oi2(CH3COO)i6(H 2 0)4] in the ground 
state (a) and the spin- flip state (b). Mnpiip indicates the 
ion whose spin direction was flipped when constructing 
the SF state. The lower panels display charge density isosur- 
faces (local density of states) respectively for the GS HOMO 
(GS1), and for both the HOMO (SF1) and the HOMO-4 
(SF2) of the SF configuration. Note that SF2 has opposite 
spin direction with respect to the remaining HOMOs around 
E F . Color code: Blue=Mn, Red=0, Green=C, light blue=H. 



We construct a Mni2-based two terminal molecular de- 
vice by sandwiching Mn X 2 in between two Au (111) sur- 
faces. The molecule is anchored to Au by a thiol group 
attached to the CeH 4 ligand, as in one of the published 
transport experiments [2 . In order to limit the lateral 
dimensions of the unit cell for the transport calculations, 
we remove 12 of the possible 16 ligands and passivate the 
remaining bonds. The final simulation cell contains 620 
atoms of which 480 are Au atoms of the electrodes, with 
the basis set expanding over 2672 local orbitals for each 
spin direction. We have carefully checked that the first 



few HOMOs and LUMOs around E^ are not affected by 
the presence of the ligands or by the passivation. Since a 
full relaxation of the entire unit cell from DFT is numer- 
ically too demanding we first optimize the geometry of 
the molecule plus ligands and then we calculate the equi- 
librium bonding distance between the ligands and the 
hollow site of the Au (111) surfaces. The resulting cell is 
presented in the inset of figure [7] (see also supplementary 
material) . 

Transport calculations are performed with the 
Smeagol code [TJ [2], which combines the non- 
equilibrium Green's function method (NEGF) [T7J [18], 
with DFT as implemented in Siesta [3]. In this case the 
rather large size of the simulation cell requires a number 
of technical solutions to boost the convergence. We have 
used our recently developed singularity removal proce- 
dure to evaluate accurately the leads self-energies [19], 
and implemented a mesh refinement scheme for integrat- 
ing efficiently the non-equilibrium charge density over en- 
ergy (see supplementary material). The entire package 
was parallelized and optimized to run on a large scale 
computational facility. 

The calculated I-V curves, their spin decomposition 
and the associated differential conductances, G(V), are 
presented in figure [7] for Mni2 both in the GS and the 
SF configuration. The most notable feature is the rather 



60 
40 



20 



% o 



-20 



-40 



-60 



1 

(a) 






— GS 




— SF ■ 


M§Ls 









-0.4 




-0.2 0.2 0.4 -0.2 0.2 

V(Volt) V(Volt) 

FIG. 2: Transport properties of a Mni2 two-probe device. 
The I-V curves for both the GS and the SF configuration 
are presented in (a), while (b) and (c) show the I-V spin- 
decomposition and the differential conductance respectively. 
The inset in panel (a) illustrates the simulation cell. 

stark difference between the GS and SF curves, which 
demonstrates that a single current readout enables one 
to distinguish between the different spin states. Impor- 
tantly this does not require the readout of the two spin- 
components of the current, i.e. a spin-polarized measure- 
ment, but it is simply obtained by comparing the general 
shape of the I-V characteristics. 

In particular the distinctive fingerprint of the two mag- 
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netic configurations is in the presence of NDRs at differ- 
ent bias positions for the two configurations. These orig- 
inate from the dynamics of the molecular energy levels 
under bias [20], their charging and their electrical po- 
larization. The main concept is that the different spin 
states of the molecule are associated to single particle 
levels with different orbital symmetry and hence with dif- 
ferent response to an external electrical potential. Thus 
a particular spin configuration translates into a distinct 
orbital configuration, which in turn is electrically read- 
able. The details of how the various molecular orbitals 
respond to the external potential are then determined by 
their coupling to the electrodes and their charging energy. 

These are investigated in figure |8j where we present the 
transmission coefficient as a function of energy, T{E), cal- 
culated at different voltages for both the GS and the SF 
configuration. In order to understand our results let us 
recall the main features of electron transport through a 
generic molecular orbital ip a . When ip a is weakly coupled 
to the electrodes the associated DOS, D a (E), is approx- 
imated [21] by 
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where j a = 7^ + 7^ and 7^ /ft (7^ / ft) is the transmis- 
sion rate for electron hopping from the left- (right-) hand 
side electrode to i/j a . The single particle energy e a de- 
pends on the charging state of the molecule. In mean 
field approximation s a scales linearly [4] with the orbital 
occupation, n a: which in turn is determined by 

J — 00 lot 

Here f L (f R ) is the Fermi function evaluated at E — /xl 
(E — /xr) with /iL (/ir) the chemical potential of the left- 
(right-) hand side electrode (/xl — Mr = e ^0- Finally the 
transmission coefficient is T a (E) = 27rj^j R D a (E) / j a . 

The expression ([7| establishes that the steady state 
current is a balance of electron fluxes to and from both 
the electrodes. Thus if a molecular level is more strongly 
coupled to one of the electrodes, as in an STM geome- 
try, its single particle energy will be pinned to the chem- 
ical potential of that electrode, and the corresponding 
transmission peak T(e a ) will move in energy as a func- 
tion of bias accordingly. In contrast, when 7^ 
then T a (e a ) as a function of V is only determined by the 
level charging energy. Crucially in both the situations 
the level broadening j a does not depend on bias and the 
width of the transmission peak does not change with V. 
As a consequence /(V) is solely determined by the posi- 
tion of the transmission peaks. Since upon increasing the 
bias charging alone cannot expel from the bias window a 
molecular level which is already inside at a lower bias, no 
NDRs are expected. This is the most typical situation in 
molecular junctions. 

In contrast the picture presented in figure [8] is rather 
different. Let us take the GS as an example. At V = 
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the first four HOMOs are less then 50 meV from E F 
and approximately 300 meV away from the first of the 
LUMO, hence they determine the low bias current. At 
positive bias the current is initially determined by the 
first HOMO, which is pinned at E F . The HOMO moves 
into the bias window at about 100 mV and then roughly 
maintains its energy position at any larger positive V. 
Notably for V < 100 mV the height of T(euoMo) is near 
to 1, suggesting that such a state is equally coupled to 
the electrodes. However when the level enters the bias 
window T(shomo) reduces dramatically (note that the 
scale of Fig. [8] is logarithmic) , indicating that now the 
coupling to the electrodes has changed. In general the 
magnitude of T(e a ) reduces if 7^ 7^ 7« • We then con- 
clude that for the HOMO the coupling to the electrodes 
becomes asymmetric with bias. 

We can extract a quantitative parameter describing the 
dynamic coupling between the electrodes and the molec- 
ular orbitals as a function of bias from the spectral rep- 
resentation of the non-equlibrium Green's function. For 
each molecular level one can evaluate the contribu- 
tions to the imaginary part of the self-energy originat- 
ing respectively from the electronic coupling to the left 
(7^) and right (7^) electrode. In addition to the total 
broadening 7^ we now can define the coupling asymme- 
try parameter r a = (7 J - 7^ + 7^ ). This is +1 
(-1) for molecular orbitals more strongly coupled to the 
left (right) contact and for coupling equally strong to 
both the electrodes. 

In figures [4] and [5] we trace e a from the position of 
the various peaks in the transmission coefficients as a 
function of V. The colors encode the value of r a (in 
Fig. [4| and j a (in Fig. [5|, and the two straight black 
lines mark the bias window. Let us start again from the 
GS [panels (a)]. As mentioned before for small V the 
first HOMO is pinned to the lower bound of the bias 
window and, confirming our initial guess, it is roughly 
equally coupled to the electrodes [Fig.|4^a)]. However for 
V ~ 100 mV it suddenly jumps into the bias window. In 
doing so a remarkable effect takes place: its coupling to 
the leads becomes extremely asymmetric, r a ~ 1, while 
7 a remains approximately constant. Considering the fact 
that the contribution to the current from a transmission 
peak, I a , is approximately 
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we conclude that as soon as the peak enters the bias 
window the total current actually decreases. This creates 
the NDR at 100 mV shown in figure [7jc). Similarly the 
second HOMO enters the bias window at 200 mV, again 
with a drastic change of r a . This produces the second 
NDR at 200 mV. 

For V < the situation is analogous although quanti- 
tatively different. Again at low bias the HOMO follows 
the lower boundary of the bias window with an almost 
equal coupling to the two electrodes. Then at about 
—200 mV it suddenly enters the bias window together 
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with the HOMO-2. Similarly to the V > case entering 
into the bias window is accompanied by a drastic change 
in coupling. This time however r a ~ — 1, indicating that 
the molecular level is now more strongly coupled to the 
right-hand side electrode. As a result the NDRs at neg- 
ative V appears only at —200 mV. 

We now turn our attention to the SF configuration. 
The most notable feature is represented by the minor- 
ity HOMO, for which the coupling to the leads is almost 
insensitive to the bias and it is rather equally strong to 
both the electrodes [see Figj4^c) and Fig[5|c)]. Therefore 
the level always follows the lower boundary of the bias 
window and remains occupied at any voltage, thus re- 
sulting in an almost linear contribution to the I-V. The 
majority component of the current instead displays a be- 
haviour qualitatively similar to that of the GS. This time 
however the first two HOMOs enter the bias window to- 
gether but only at negative voltages (around —100 mV), 
again in correspondence to an increased asymmetry of 
the electronic coupling with the electrodes [r a ~ —1, see 
Fig.gb)]. This generates the NDR at -100 mV. 

Finally we answer the following question: what is the 
mechanism behind the bias-dependence of r a ? Let us 
consider the isosurfaces of Fig. |8j displaying the elec- 
tronic wave function of a given molecular level in pres- 
ence of V. For example let us look at the panels 2, 4 
and 6 for the GS. These represent the wave- function of 
the first HOMO respectively at 300 mV, and -300 mV. 
Clearly as bias is applied, the wave function experiences a 
substantial polarization, with the electron clouds chang- 
ing from a uniform distribution at V = to one which is 
considerably localized either to the left- or to the right- 
hand side of the junction. The molecular orbital dis- 
tortion under bias originates from molecular orbital re- 
hybridization and always involves a number of molecular 
states. For instance, when the HOMO and HOMO-1 
enter the bias window at V > the re-hybirdization in- 



volves the first several HOMOs, as one can deduce from 
the fact that for most of them there is a change in r a [see 
Fig. [4ja)]. Importantly for the SF configuration there is 
only one minority spin level around Ep, this cannot hy- 
bridize with any other state (unless a spin-flip event oc- 
curs), and therefore its coupling to the electrodes changes 
little with bias. This fundamental aspect is what creates 
profoundly different I-V curves in the two cases. 

Crucially orbital re-hybridization results in the molec- 
ular levels entering the bias window always to display 
largely asymmetric coupling to the electrodes. This al- 
lows them to enter the bias window and to conduct with- 
out charging the molecule, which remains neutral at any 
bias investigated. We therefore conclude that the ob- 
served NDRs are features in the I-V occurring at energies 
lower than the typical charging energy of the molecule. 
Note that in general orbital re-hybridization may not 
only lead to NDRs, but also to a sharp increase of G(V). 
Importantly these are always low-energy satellite features 
to the main Coulomb blockade signal. 

In conclusion we have demonstrated that the I-V char- 
acteristics of Mni2-based devices are sensitive to the in- 
ternal magnetic state of the molecule, so that an electrical 
readout can distinguish between two different magnetic 
configurations without the need of resolving the two spin 
components of the current. Such a sensitivity originates 
from the bias dependence of the electronic coupling be- 
tween the molecular levels and the electrodes. In the 
case examined here this provides a novel mechanism for 
NDRs, which indeed are the final fingerprints of a spe- 
cific magnetic state of the molecule. Significantly this 
represents also a new mechanism for NDR, requiring a 
number of closely spaced molecular levels close to the 
electrodes' Fermi level. Single molecules magnet, with 
their d-electron manifold, appear as the ideal material 
system for such an effect. 
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FIG. 3: Transmission coefficient, T(E) as a function of bias for both the GS (left) and the SF configuration (right). The 
isosurfaces display the wave- functions of a number of molecular levels at different V. For these we also report the value of both 
r a and j a . Note the dramatic re- hybridization with bias leading to highly asymmetric wave-functions. 
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FIG. 4: Trace of r a as a function of energy and bias. The colored lines denote the position of the peaks in the transmission 
coefficient and the color encodes the magnitude of r a . The two solid straight black lines mark the boundaries of the bias 
window. Panel (a) is for the GS and majority spins, panel (b) is for the SF configuration and majority spins, panel (c) is for 
the SF configuration and minority spins. Color code: red r a = 1, green r a = 0, blue r a — — 1. 
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FIG. 5: Trace of j a as a function of energy and bias. The colored lines denote the position of the peaks in the transmission 
coefficient and the color encodes the magnitude of j a . The two solid straight black lines mark the boundaries of the bias 
window. Panel (a) is for the GS and majority spins, panel (b) is for the SF configuration and majority spins, panel (c) is for 
the SF configuration and minority spins. Color code: red j a = 10 -3 eV, blue j a = 6.5 • 10 -5 eV. 
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Supplementary Materials 
I. COMPUTATIONAL METHOD 




FIG. 6: Orbital resolved density of states around the Fermi 
level (purple vertical line) for [Mni 2 Oi2(CH3COO)i6(H 2 0)4] 
in the S = 10 ground state. Note that the HOMO-LUMO 
gap is an intra-d gap with states originating mainly from Mn 
ions (mainly 3+, see main text). Note also the strong Mn-0 
hybridization both in the HOMO and LUMO. 

Calculations have been all performed with the DFT and 
NEGF smeagol code [IJ [2]. We use standard scalar rela- 
tivistic pseudopotentials with the following reference config- 
urations: H Is 1 , C 2s 2 2p 2 , S 3s 2 3p 4 , O 2s 2 2p 4 , Mn 4s 2 3d 5 
and Au 6s 1 . The atomic basis set is constructed as follows 
C: DZ-s, DZ-p; H: DZP-s; O: DZ-s, DZP-p, SZ-d; S: DZ-s, 
DZP-p, SZ-d; Mn: DZP-s, SZP-p, DZ-d; Au: DZ-s (the no- 
tation is SZ=single zeta, DZ=double zeta, P=polarized)[3]- 
The real space grid has an equivalent cutoff energy of 400 Ry. 
The charge density is obtained by splitting the integral of the 
Green's function (GF) into a contribution calculated over the 
complex energy plane and one along the real axis 1 . The 
complex integral is performed over a uniform mesh of 512 
imaginary energies, while for the real part we have imple- 
mented a new mesh refinement algorithm, necessary to in- 
tegrate the extremely sharp features of the DOS. Such an 
algorithm consists in evaluating the integral over an initial 
coarse energy mesh after the GF has been artificially broad- 



ened. Then peaks in the GF are detected and a denser energy 
mesh is generated around them. The next approximation is 
then calculated over this new grid after the artificial broad- 
ening has been decreased. Such a procedure is repeated un- 
til no artificial broadening is left. Typically the final mesh 
includes 1000 energy points with the denser energy spacing 
being around 10 -5 eV. 

The hopping rates of a molecular orbital to the electrodes, 
7^ and 7^, used in the analysis of the transmission coeffi- 
cients, are calculated by using the spectral representation of 
the GF associated to the scattering region 



(4) 



where ip a and ip a are the right and left eigenvectors of if e ff 
and {soi — are the associated eigenvalues. Finally 



(5) 



with Hm the Hamiltonian of the scattering region (molecules 
plus part of the leads) and Sr (Sl) is the self-energy for the 
right- (left-) hand side electrode pQ. The coupling constants 
la ( n — L 5 R) are obtained as 



l2 fi[^n - S^)]^ a 



(6) 



One can easily demonstrate that with this definition j a = 
+ la an d the orbital occupation of is [equation (2) of 
the main text]. 
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FIG. 7: Orbital resolved density of states around the Fermi 
level (purple vertical line) for [Mni 2 Oi2(CH3COO)i6(H 2 0)4] 
in the S = 10 ground state. Note that the HOMO-LUMO gap 
is an intra-d gap with states originating mainly from Mn 3+ . 
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II. MNi 2 DENSITY OF STATES 

In figure [6] we present the orbital resolved DOS for 
[Mni 2 Oi2(CH 3 COO)i6(H20)4] as calculated with DFT-GGA 
for the S = 10 ground state. We note that the HOMO-LUMO 
gap is an intra-d gap formed mainly between molecular or- 
bit als with amplitude over the Mn ions. These have a large 
degree of hybridization with O and they are fully polarized 
around the Fermi level (Ef). Also note that there is little 
contribution to the DOS around Ef from the C p orbit als. 

In figure 17] we further analyze the DOS by projecting it 
over the Mn" and Mn 4+ ions individually. The figure demon- 
strate that the HOMO-LUMO gap is indeed solely determined 
by Mn 3+ d-shell with the first of the Mn 4+ states being at 
least 1 eV away from Ef. 



III. TRANSPORT SIMULATION CELL 




FIG. 8: Transport simulation cell used in this work. Color 
code: Blue=Mn, Red=0, Green=C, Light Blue=H, Yel- 
low=Au, Dark Yellow=S, Purple=Mn (nipped). 



The transport simulation cell is constructed from a 
[Mni20i2(CH 3 COO)i6(H 2 0)4] molecule comprising 16 thiol- 
terminated C6H4 ligands. We cut and passivate 12 of those 
16 ligands in order to reduce the lateral size of the cell. Such 
final molecule is relaxed in vacuum by standard conjugate 
gradient method and then positioned at a minimal energy lo- 
cation on the Au (111) surface. The search for the minimum 
is conducted as follows. First we place the S atom of one of 
the remaining 4 ligands at the hollow site of the Au (111) 
surface [this is one of the preferential bonding sides for thiol 
on Au (111) 01] j and then we rotate the molecule about this 
position in the search of the total energy minimum. The fi- 
nal device simulation cell (see figure Fig. [7| also includes five 
Au atomic planes on each side of the molecule. The lateral 
dimensions of these are those of a 6 x 4 supercell constructed 
from the primitive cell along the fee (111) direction. The unit 
cell for the electrodes is also a 6 x 4 supercell. The resulting 
scattering region (where the self-consistent calculation takes 
place) is presented in figure u\ Note that two of the 12 Mn 
atoms are colored in purple (the other ones are blue). These 
are the ions whose spins are flipped in forming the spin-flip 
configuration of the molecule. 
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